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Abstract 

The nonlinear interactions that evolve between a planar or nearly planar Tollmien- 
Schlichting (TS) wave and the associated longitudinal vortices are considered theoreti- 
cally, for a boundary layer at high Reynolds numbers. The vortex flow is either induced 
by the TS nonlinear forcing or is input upstream, and similarly for the nonlinear wave 
development. Three major kinds of nonlinear spatial evolution, Types I-III, are found. 
Each can start from secondary instability and then becomes nonlinear, Type I prov- 
ing to be relatively benign but able to act as a pre-cursor to the Types II, III which 
turn out to be very powerful nonlinear interactions. Type II involves faster stream- 
wise dependence and leads to a finite-distance blow-up in the amplitudes, which then 
triggers the full nonlinear three-dimensional triple-deck response, thus entirely altering 
the mean-flow profile locally. In contrast, Type III involves slower streamwise depen- 
dence but a faster spanwise response, with a small TS amplitude thereby causing an 
enhanced vortex effect which, again, is substantial enough to entirely alter the mean- 
flow profile, on a more global scale. Streak-like formations in which there is localized 
concentration of streamwise vorticity and/or wave amplitude can appear, and certain 
of the nonlinear features also suggest by-pass processes for transition and significant 
changes in the flow structure downstream. The powerful nonlinear 3D interactions II, 
III are potentially very relevant to experimental findings in transition. 
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1. INTRODUCTION 


Experiments in boundary-layer transition tend to show that, depending on the in- 
put amplitude and spectrum of the disturbance upstream, a two-dimensional Tollmien- 
Schlichting (2D TS) incident wave readily succumbs to 3D secondary instabilities down- 
stream which then induce fully 3D nonlinear responses in the subsequent transitional 
boundary-layer flow: see, e.g., Arnal’s (1984) review. The resonant-triad nonlinear mech- 
anism (Craik 1971, Smith and Stewart 1987) provides one theoretical explanation for 
this nonlinear 3D interaction for low input amplitudes. Also on the theoretical front, for 
boundary-layer transition, is the related work by Hall and Smith (1984, 1988a) on oblique 
TS input or incident vortex flow with nonlinear interaction, while Benney and Lin (1960) 
make other interesting suggestions and Herbert (1984) presents interesting accounts of 
the first, but only linear, secondary 3D phase: see also the criticism in Hall and Smith 
(1988a) and later comments. Again, there are recent analyses tackling the corresponding 
channel- flow nonlinear interactions by Srivastava and Dallmann (1987), Hall and Smith 
(1987, 1988b, c). Our aim in the present theoretical study is to derive the scales, structure, 
governing equations and solution properties for nonlinear 3D interactions set up between 
incident 2D, or nearly 2D, TS waves and their induced vortex motion, in incompressible 
boundary-layer transition. 

As TS waves are involved it is possible, and in some ways natural, to start with 3D 
triple-deck theory (we extend this subsequently) since that is known to identify in a rational 
way the framework of 2D and 3D linear and nonlinear TS waves (Smith 1979, Hall and 
Smith 1984) near the first, lower-branch, neutral station, say for the Blasius boundary layer 
on a flat plate. This theory however also incorporates directly the induced longitudinal- 
vortex motion provided the latter’s streamwise length scale is not excessive, as described 
subsequently. The three regions of the triple-deck structure are the lower, main, and upper 
decks. The lower deck lies closest to the plate surface, inside the boundary layer, and the 
flow there responds to a 3D nonlinear unsteady-inertial-pressure gradient-viscous force 
balance. This deck coincides with the critical layer for very small, linearized, disturbances, 
although our emphasis here is necessarily on the nonlinear range. Above the lower deck the 
motion in the main deck spanning the boundary layer is displaced simply in a quasi-steady 
planar fashion, due to the lower deck, and the displacement effect is thereby transmitted to 
the upper deck lying outside the boundary layer. There the inviscidly provoked pressure 
response has to be consistent with the wall pressure driving the lower-deck flow, thus 
producing viscous-inviscid interaction. The whole nonlinear evolution is governed mainly 
by the flow features in the lower deck, where the nondimensionalized velocity components 
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and pressure have the form 


[eA*u,e 3 A*v,eA*iu,e 2 A2p] 


(l.lo) 


in terms of the nondimensionalized local coordinates and time 

[x 0 + e 3 \-*X,e B \-zY,z 0 + e a \~* Z,e 2 \~*T\. (1.16) 

Here Re = e~ s is assumed large, in tune with the large experimental values of the global 
Reynolds number Re of interest, and our attention is focussed initially around a typical 
0(1) station x = x Q > 0, z = z 0 on the plate Y = 0, x > 0, the undisturbed flow being 
predominantly in the x direction and giving an 0(1) skin-friction factor A locally at x 0 , z Q . 
With (l.la,b) holding, the Navier-Stokes equations reduce to the 3D unsteady interactive 
boundary-layer equations for u, v, w, p as functions of X, Y, Z, T in the lower deck, subject 
among other conditions to matching with the displaced main-deck solution. The unknown 
displacement oc —A involved there is related to the unknown pressure p in (1.1a) by the 
upper-deck properties applying outside the original boundary layer. 

This nonlinear interactive system is addressed in Section 2 with regard to vortex/TS 
interaction at reduced amplitudes; three principal types (I-III) of interaction are identi- 
fied, in fact, and we choose to concentrate first on “Type I.” If the input amplitudes are 
sufficiently low, on the other hand, nonparallel-flow effects due to the slow variation of 
A with x matter considerably and structural scales different from those in (1.1a, b) come 
into play, as described in Section 3. These scales are derived from an adaptation of those 
above; alternative derivations come from the related Hall and Smith studies (1984, 1988a, 
and work in progress). Another point of note here is that our concern is with spatial 
evolution. This seems to tie in better with the experimental and real-life behavior than 
does a temporal-instability analogy of the sort used by Herbert (1984) for example, who 
also makes the irrational shape-factor assumption for the 2D TS contribution as opposed 
to the present rational approach based on nonlinear evolution equations for both the TS 
wave and the induced vortex pattern. 

The flow properties resulting from the vortex/TS nonlinear interaction equations de- 
rived in Sections 2,3 are considered in Section 4. The flow-structural analysis and results 
in Sections 2-4 are specifically for “Type I” interactions in which the typical streamwise 
(X) variation is relatively slow, as is the spanwise (Z) variation associated with an almost 
2D input upstream, producing a “warped” effect. Faster streamwise responses lead to the 
“Type II” nonlinear interaction discussed in Section 5, while faster spanwise dependence 
leads to the “Type III” interaction which is described in Section 6. Further discussion 
of these and other aspects is presented in Section 7. The Type I vortex/TS interaction, 
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for instance, which we study first, is found to provoke either a finite-distance breakdown 
or a long-distance sustained effect, both of these producing a marked alteration in scales 
and hence structural change downstream and possibly leading on subsequently to Type II 
and Type III interactions. Further, the Type II, III interactions turn out to be much more 
powerful than the Type I interaction and suggest (e.g.) streak-like formations downstream, 
associated with local spanwise concentration of the streamwise vorticity and the TS am- 
plitudes. Both Types II, III are potentially very relevant to the experimental findings as 
(unlike other theories) they lead to complete alteration of the mean-flow profiles. 


2. THE TRIPLE-DECK VERSION FOR WARPED TS-VORTEX INTER- 
ACTION (TYPE I) 


We start by posing the triple-deck problem which governs the linear or nonlinear evo- 
lution of 2D and 3D TS waves and their mean-flow effects, initially at least, and which 
stems from substitution of (l.la,b) into the Navier-Stokes equations. This requires us to 
tackle the unsteady 3D nonlinear interactive boundary-layer equations 

u x + v Y + wz = 0, (2.1a) 

< Ut + UUx + VUy + WUz — — Px(X, Z, T) + Uyy, (2.16) 

, Wt + UWx + VWy + WWz = ~Pz{X, Z , T) + VJyy, (2.1c) 


describing the lower-deck response, in scaled form, for the unknowns (u, v,w)(X, Y, Z,T ), 
(p, A) (X, Z, T) subject to the boundary conditions 


u = u = u> = 0at Y = 0, 
j u ~ Y + A{X, Z, T ) , w — ► 0, as Y — ► oo, 

°° r a^A({^,T)/ded(d n 




[( X - W + (z - r,)>}i 


(2. Id) 
(2-le) 

( 2 - 1 /) 


Here the normal momentum balance reduces to the usual requirement that dp/dY — 
0 as assumed in the above, and (2.1d-f) represent in turn the no-slip condition at the 
wall, the displacement effect (a —A, unknown) on the main-deck motion further from 
the wall, and the pressure-displacement interaction with the flow outside the boundary 
layer, produced via the upper deck where linearized potential-flow theory holds, yielding 
the double Cauchy-Hilbert integral relation shown. This relation can be replaced by the 
complete upper-deck formulation of solving Laplace’s equation for the pressure, p say, 
with p zero in the farfield (X 2 + y 2 -*• oo) and equal to p(X, Z,T) at the wall (y 0+), 
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where also the gradient dp/dy is to equal d 2 A/dX 2 {X, Z, T ), in terms of the upper-deck’s 
scaled normal coordinate y. The nature of the triple-deck flow in (2.1a-f) appropriate to 
nonlinear interaction between a warped 2D TS input wave and its induced vortex flow is 
discussed here, with the next section then considering the required extension of the theory 
to encompass global nonparallelism effects. 

The starting point is that the interactive triple-deck system (2.1) captures traditional 
linear TS waves (Smith 1979) as small unsteady perturbations of the basic boundary-layer 
shear flow u = Y,v = w = p = A = 0. The waves are stable or unstable for a prescribed 
scaled input frequency 0 according as n^flc(« 2.30), where the critical frequency fi = fl c 
defines the lower branch of the neutral curve for the boundary layer. Near the critical 
frequency, or equivalently, near the critical station z 0 = x c for an imposed dimensional 
frequency, interplay usually occurs between nonlinear amplitude-dependent effects and the 
relatively weak initial growth/decay present, for a small 2D or 3D disturbance. Supposing 
the specified input disturbance amplitude to be h, relative to the normalized form (2.1), 
and the input wave to be almost 2D, with a small “warping” factor (3 say, so that the 
characteristic variation has d z ~ /?, we argue by orders of magnitude as follows (see also 
Fig. 1). The vortex flow produced by the small warped-TS wave arises from inertial forces 
as an amplitude-squared response but with the tu-component of the vortex multiplied also 
by /?, since the typical TS tu-velocity is reduced to 0 (h/3), thus reducing the the inertial 
secondary forcing to 0(h 2 /3) via contributions such as udw/dX. That 0 (h 2 0) u/-component 
in the vortex flow grows logarithmically with large Y, however, due to the decreased inertial 
reaction of the vortex. So an outer buffer zone is brought into play. This has F-extent 
of order k~* to bring about an inertial-viscous balance for the vortex, i.e. balancing Yk 
against dy because of the basic shear (u « Y) and the unknown slower X-variation dx ~ k 
(say) in the vortex flow due to the slight warping present. In the buffer zone the vortex 
velocity w is still 0 [h 2 0), apart from a logarithmic factor, and so continuity suggests a 
vortex u-component of size (3 2 h 2 /k (giving ux ~ wz ~ (3 2 h 2 ) and hence a shear-correction 
effect du/dY of order (3 2 h 2 /ki. Our reasoning then is that a sensitive interaction is likely 
to take place when the three relative corrections, 0(/? 2 ) from the warping, 0(/i 2 ) in the 
traditional nonlinear amplitude-cubed feedback, and 0(f3 2 h 2 /&»), due to the nonlinearly 
induced vortex shear, are comparable. So the regime 

/?~/i,A;~h 3 (2.2) 

is indicated as a central one, which we take as defining “Type I” interactions. Other 
regimes of the warping factor (3 and the streamwise-variation factor k can be examined 
but they may be regarded either as extreme limits of that in (2.2) or as distinct ones. In 
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fact a matter of significance at this stage concerns the small order h s in (2.2) corresponding 
to the Type-I streamwise variation. This is formally less than the traditional influence of 
relative size h 2 usually associated with amplitude-cubed nonlinear spatial evolution in a 
purely 2D input wave. That aspect leads on to the study of the alternative Type II and 
Type III interactions as described subsequently in Sections 5,6. 

(a) The near-wall zone . 

Guided by the ordering arguments above, the pattern of the flow solution in the near- 


wall region where Y is 0(1) has 

u = Y + hUi + h 2 {U 2 + A 9 Y) + h s U z + • • • + h 4 U 4 + h 5 U 5 + •■■, (2.3 a) 

v = hV 1 + h 2 V 2 + h s V s + h 4 V 4 + h s V s + • • • , (2.36) 

w = h 2 Wi + h s W 2 + h*W z + h*W 4 + h e W 5 + ■ • • , (2.3c) 

p = hPi. + h 2 P 2 + h 3 P 3 + h*P 4 +■■■, (2.3 d) 

A — hAi h 2 A 2 + • • ■ , (2.3e) 

along with the extra scaling 

dx &x + ^ S ^x + ' “ hdj. (2.4) 


Here ( U 2 , Vj, W 2 , Pi) represent the slightly warped 2D TS input wave, (U 2 , V 2 , P 2 ) are mean- 
flow correction and second-harmonic effects that are quasi-planar, similarly (I/3, V 3 ,P 3 ) 
include the amplitude-cubed forced TS response, while the induced-vortex motion is rep- 
resented in the terms A 3 y and U 5 , V 6 ,W 2 and is independent of the faster scales X, T. The 
induced shear A 3 (X, ~Z) is an unknown feedback effect from the outer buffer zone which is 
studied subsequently in (b). In addition, a streamwise scale X\ = h 2 X for the modulation 
of the TS amplitudes due to essentially 2D nonlinear forcing is absent here (except in the 
phase, see below), leaving the slower scale X(= h s X) to operate on the TS modulation 
in which the nonlinear TS /induced-vortex interplay is set up. An assumption of near- 
criticality, e.g. D = n c + 0 (h 2 ) or Xo = x e + 0(h 2 ) in unsealed terms, is also implicit here in 
view of the relatively small amplitude growth present initially. The disparity in certain of 
the velocity scales in (2.3), e.g. in the main TS contributions Ui, Vi,Wi and the induced 
vortex contributions U$,V 3 ,W 2 , is caused by the input warping defined in (2.2), (2.4). 

The dominant governing equations that we need to address are obtained formally from 
substitution of (2.3), (2.4) into (2.1) and are the following: first, for continuity 

Uix + V\ y = 0, (2.5a) 
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U 2 x + Yiy — 0, 
Uzx + V 3 y + = 0 ; 


(2.56) 

(2.5c) 


second, from the X-momentum balance, 


U\ t + YUxx + V X = —Pix + Uiyy, 


(2.6a) 


U 2T + YU 2X + [tWuc + VJJiy] + V t = -P 2X + U 2Y y, (2.66) 


UzT+YUsx+^YUix+lUiUix+UzUix+ViUzY+ViUiY^Vs+ViXs+WiUtf = — Psx+U 3 yy] 

(2.6c) 


and, third, in the ^-momentum balance, 


W 1t + YW 1x = -P 1 j+W 1Y y, 


(2.7 a) 


[UiW lx + VMy] = W 2YY . (2.76) 

In (2.7b) , only the induced-vortex terms, which are independent of the fastest scales X , T, 
are retained. Further, the successive governing equations in (2.5) - (2.7) show the grad- 
ual intrusion of three of the four major influences present, the warping through the ~Z 
variations, the amplitude-squared, -cubed, -etc., nonlinear responses and the feedback 
induced-vortex-shear contributions through the A s(X,Z) factors, but not yet the stream- 
wise modulation in X which acts more significantly in (b) below. 

The velocity and pressure fields above can be written now in component form, 


U\ — U n E + c.c, U 2 — {U 22 E 2 + U 2i E) + c.c. + U 2 Q t 
Uz = U 33 E 3 + -",U 4 = U«E 4 + -‘-, 


( 2 . 8 ) 


and so on, and the solutions for the various components are obtained in sequence. Here 
the powers of 

E = exp(t(aX - flT)] (2.9) 

present contain all the fast dependence, with the main wavenumber a and frequency fi 
being assumed real (see (2.12) below), so that the unsteady flow solution is near neutral. 
Also c.c. or, later, an asterisk denotes the complex conjugate function. The dominant, 
almost planar, TS wave is then controlled by the equations and constraints 


iaU n + V 11Y = 0, 

— iClUu + iaYUn + Vji = —iaPu + t/nyy, 
U\i = Vn = 0 at Y = 0, U\\ — ♦ An Y — ► 00 , 
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(2.10a) 

(2.106) 
(2.10c, d) 


P n = ccA n , 


(2.10e) 


from (2.5a), (2.6a), with (2.1d-f), and (2.10a-d) yield the shear-stress solution 

U llY = B u Ai(€),t = (*a)*(y - n/a), 

► 

Bu = (ta)fP n /Ai' 0 , (ia)iP U K = A n Ai' 0 


( 2 . 11 ) 


where the subscript zero refers to evaluation at £ = &> = — (ta) $f]/ a, i n denotes exp(m7r/2) , At 
is the Airy function, and k = Ai(s)ds. Hence the pressure-displacement law (2.10e) 
leads to the eigenrelation 

t’s’a? = Ai'q/k (2.12a) 

between a, fi. Both of these are real as required only for the values £o = —diis,Ai' 0 /K = 
d 2 i*\di « 2.30 , ^2 « 1.00], i.e. the values 

a = df,fi = djdf (2.126) 


are fixed. The fundamental pressure amplitude Pu(X, Z, • • •) remains undetermined at this 
level: see also (2.14) ff below. In addition, the spanwise velocity component Wn associated 
with the incident TS warping is only a passive response so far, given by 

Wu = P ll i£(()/(ic')i (2.13) 


from (2.7a), where £ is the solution of £" — ££ = 1, Co = -C(oo) = 0. 

Moving on to the second-order responses, we may proceed similarly to determine 
^ 20 , 21,22 from the component equations inferred from (2.5b), (2.6b), which again are quasi- 
2D. The working for the forced second-harmonic TS term U 22 and the forced mean-flow 
correction U i0 has been done before, however, by Smith (1979), giving in particular V 2 o> P 20 
identically zero but U 2 0 , A 2 o nonzero. The un-forced extra fundamental which may be pro- 
duced here for compatibility at higher order simply has (U 2 i,V 2 i, An) equal to P 21 /P 11 
times (Uii,Vn,Au) with P 21 still unknown. 

The third-order working for l7s,Vs components is then as in Smith (1979) except that 
they feel the additional influences of the input warping, in the two d/dZ terms in (2.5c), 
(2.6c), and of the induced-vortex flow, in the two A 3 terms in (2.6c). With the warping and 
vortex influences incorporated, then, the P-components contained in those two equations 
lead, through a compatibility requirement, to an amplitude-modulation equation part of 
which can be picked out from the above paper or from Hall and Smith (1984). We find 
therefore the modulation equation 

iaP n + cA 3 Ph + d-=± = eP n + Mi|Pn| 2 (2.14) 

dZ 
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for the unknown TS pressure amplitude P u , on use of (2.13). The coefficients involved 
here, 

c = 2Di£oa 2 Aio/3A* Ai' 0 — 5/3, d = 1/8 a 2 , (2.15a, 6) 

e = ia 1 a, f = a x a, a = — t(c + 3)/a, (2.15c, d,e) 

are in keeping with those in the previous studies mentioned earlier, for small /?, with 
(complex) representing the effective wavenumber shift and a x the Stuart-Landau constant 
for the pressure. Also D = 1 + /cfo/Atg, A = ia. Clearly, the effects of the induced-vortex 
motion and the warping on the TS wave are represented by the terms in A 3 and d 2 jdZ in 
(2.14). The other term on the left in (2.14) is the relatively fast phase effect noted earlier, 
corresponding to an extra exp(iaXi) factor necessary in Pn where the phase a is real and 
constant, or associated with a contribution h 2 ia added to the streamwise variation in (2.4). 
This does not affect directly the vortex properties below and indeed <7 could be absorbed 
into the term e in (2.14). 

Lastly here, the driven vortex motion in (2.7b) needs to be addressed. The relevant 
P°-components give the equation 

W 20YY = ia{-U u w ; 1 + UiiWu) + V n w; lY + V^W llY (2.16 a) 

for W20, which is subject to 

W 2 o = 0 at Y — 0, W 2 oy — * 0 as Y — ► 00 . (2.166, c) 

Due to the forcing terms in (2.16a) decaying algebraically like Y~ 2 at large Y (from (2.13)), 
the integration for is found to produce the logarithmic-growth result 

W 20 - d T {\ P^Dl-lnY + (f>]{as F — ► 00 ) (2.16d) 

where the constant <t> is to be determined numerically from (2.16a-c). The behavior (2.16d) 
and the corresponding growth a YtnY,Y 3 inY in V 50 , U$o are responsible for the strength 
of the induced-vortex flow set up in the outer zone discussed next. 

(b) The outer-buffer zone. 

In the buffer zone Y is larger, Y = h~*Y with Y of order unity, and the flow solution 
is of the form 

u = h 1 F + h(Ai + U 3 ) + • • • + h 2 A 2 + • • • + h 5 RiY + • • • (2.17a) 

v — — A\xY + h[Qi — A 2 xY) + • • • + h 3 v 3 + • • • , (2.176) 
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w = h"* {tnK)w 3£, + h?{C]Y + UJ3) + • • • , (2.17c) 

with the pressure p and displacement — A remaining as in (2.3) of course. Here the unknown 
induced-vortex contribution of interest is represented by (u s , v s ,w 3 ), dependent mainly on 
X, Y, ~Z\ although logarithmically larger vortex velocities [w Z l and similarly in u, v\ are 
also provoked, we note, these react only passively to the forcing present in the matching 
constraints for Y — ► 0+ due to (2.16d), analogous to those in the related studies mentioned 
earlier. Further, the momentum balances in (a) above show that Q 1 = —P\x — Ait and 
Ci, R\ satisfy Cix = -Pjj, R\x + C& = 0 , i.e., 

Qn = — + inA u ,iaC n = —P n j,iaRu + C n j = 0, (2-18) 

in component form. 

Most of the leading-order balances that result from substitution of (2.17) into (2.1) are 
trivially satisfied but the 2S°-components require the vortex motion to satisfy 


U 3 X + v 3 7 + W 3 J ~ °> (2.19a) 

y U 3 X "h ^3 = u $ yy> (2.196) 

^ W 3 X + ^j(l^°n| 2 )/^ = FF> (2.19c) 

after some manipulation of the forcing terms involving (2.18). The boundary conditions 
on (2.19a-c) are 

U3 = v 3 = 0 at Y = 0, (2.19d) 

u>3 ~ 9j()Pii | 2 )[— InY + 4>\ as Y — >■ 0+, (2.19e) 

U 3 Y — ♦ 0>^3 — »• 0(7 3 ) as Y — y 00, (2.19/) 


where (2.19d,e) merge the solution with the near-wall vortex form in (a), e.g. in (2.16d), 
while (2.19f) achieves the outermost displacement behavior of (2.1e) as required. In con- 
trast with the mostly viscous response of the vortex motion in zone (a), i.e. in (2.16a), the 
vortex here in zone (b) reacts in a viscous-inviscid fashion through (2.19a-c), in response 
to the warped-TS-forcing ( oc |Pn| 2 ) acting in (2.19c,e). The forced vortex flow in (2.19a-f) 
produces the (unknown) wall-shear contribution referred to earlier, 

A3 = dyu z at Y = 0+, (2.20) 

which then helps to drive the TS amplitude P n (via (2.14)) which in turn drives the 
vortex flow (2.19), and so on, thus producing nonlinear interaction between the warped 
TS disturbance and the induced vortex. 

The nonlinear TS/vortex interaction found, as represented by (2.14), (2.19a-f) [with 
(f> determined from (2.16a-d)], is studied further in Section 4, after the extension of the 
theory which is described next. 
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3. EXTENSION OF THE THEORY TO INCORPORATE NONPARALLEL- 
FLOW EFFECTS (TYPE I) 

The influence of the nonparallelism of the basic boundary layer is secondary within 
the scope of the triple-deck scales leading to (2.1) because the latter’s streamwise length 
scale [see (l.l)] is typically 0(e 3 ), much less than the 0(1) length usually associated with 
the basic nonparallelism. The difference diminishes however for the reduced 0(/i)-sized 
disturbances considered in Section 2, since their main modulation length becomes extended 
to 0 (h~ s ) times the e 3 scale, from (2.4), while the characteristic nonparallel-flow length 
scale contracts, to 0(h 2 ), in view of the shift away from the neutral position. Hence, 
formally, the two length scales can become comparable if the disturbance amplitude h is 
reduced to 

h ~ e m ,m = 7, (3.1) 

with significant nonparallelism then entering the fray. Strictly, the setting of h as an inverse 
power of the Reynolds number here requires us to start again from scratch, replacing 
(l.la,b) by a new expansion as implied by (2.3) with (3.1) inserted and similarly in the 
main and upper decks, as well as in the extra buffer zone indicated in Section 2. But 
the only new substantial effect to come into play then is that corresponding to the global 
nonparallelism (involving Ai = dX/dx at x 0 )* So we may cut matters short simply by 
adding in that effect to the previous working, along with (3.1), in retrospect. Here again 
the additional term involved may be picked out from the Hall and Smith (1984) study, so 
that the new equation for the warped-TS amplitude is 

o2 p 

ioP\\ + ( bX + cAs)Pn + d — = ePu + /Pn|Pn| 2 (3.2) 

dZ 

when (3.1) holds, with (2.19) remaining intact. In (3.2) b = Aic captures the nonparallel- 
flow effect. The TS/vortex interaction is thus controlled by (2.19), (3.2) in the regime 
defined by (3.1), as distinct from (2.19), (2.14) which apply for the higher-amplitude 
regime where 1 >> h » e m in effect. 

4. NONLINEAR INTERACTION PROPERTIES (TYPE I). 

The warped-TS/induced-vortex interaction of Type I is governed by the nonlinear 
linked p.d.e. system implied by (2.19), (3.2): 

w yy ~ Yw x = dj[\P\ 2 )/Y\ 

Tyy ~ Ytjr = ~Wj 
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(4.1a) 

(4.16) 


subject to 


(4.1c) 

(4-l<0 


w'~a 2 -(|/ , | I ){-«»y + ^},r r ^o, >sy->o, 

t — ► 0, W — * 0, as Y — ► oo, 

for the vortex, coupled with the TS amplitude P which is controlled by 

a d 2 P 

0 = (ci + c^\\X + c 3 A 3 )P + C4 2 + c 5 P|P| 2 

dZ 

where the vortex skin-friction factor A 3 is 

A 3 = r(X,0, Z), 


(4-lc) 


( 4 - 1 /) 


and with <t> determined from integration of (2.16). The coefficients involved are given by 
Ci a (e — ia) and 


C2 — c 3 — — 0.375 -f* 0.117*, C4 — — 0.0156 + 0.0465* 
c 5 = -0.719 - 0.766*, <f> = 4.15, 


(4.2) 


from a computation of (2.15), (2.16) with (2.21b). Also we have introduced t = dU/dY 
to obtain the W — t formulation above instead of the U — V — W version earlier. Here 
U,V,W are related to be velocities u 3 ,t> 3 ,u> 3 of Section 2, while P is related to Pn, and in 
the context of Section 2 the factor A x is replaced by zero, thus suppressing the nonparallel- 
flow effect, as opposed to positive or (more likely) negative values of Ai which incorporate 
nonparallelism as in Section 3. 

Computational studies of (4.1a-f) with (4.2) for various starting conditions upstream 
were made by use of a spectral method analogous to the treatment in Hall and Smith 
(1988a) except that for the vortex part (4.1a-d) we applied an implicit scheme like that in 
Hall and Smith (1988c) and the updating of P through (4.1e) was done differently, via a 
Newton iterative procedure. An alternative formulation in terms of an integral equation 
for A 3 (see Appendix A) was treated in the same iterative way, yielding similar trends. 
The results are summarized later on. 

We consider analytically now possible ultimate forms of the nonlinear interaction as X 
increases. It is useful to express P in polar form first, P = r exp(*0) with r, 9 real, so that 
(4.1e) is replaced by the two real equations 


0 — (cj r + c^AiX + c 3r A 3 )r + c^ r {r zz r ^) c 4 i{^ r 'z^ z r @zz) c 5r r > (4.3a) 

0 = (cx< + c 2 ,Ax X + c 3 ,A 3 )r + c 4t (r^ - rOj) + c 4r (2r^% + rOj^) + c 5 ,r 3 (4.36) 
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linked with (4.1a-d,f) by |P| 2 = r 2 . Then one possibility (option 1) suggested by related 
studies is that of a nonlinear exponential interaction continuing to downstream infinity, 
giving the response 

M) ~ ~ e 2, ^(lV 00 ,r 00 , A 3oo ) (4.4) 

as X — * oo, with s an unknown positive constant. From (4.3) ra,,#**, must satisfy 

0 C3 r A 3oo "b ®5r^oo> (4.5fl) 

0 = C 3 jA 3oo (4.56) 

where A 3oo oc (r^)jj is r x at Y = 0+ and r^, , W 0 0 are controlled by forced Airy equations, 
the solutions of which can be written down explicitly. The coefficients involved in (4.5a,b) 
[see (4.2)] however are such that (4.5a,b) admit no positive solutions for and 

hence the nonlinear exponential response seems to be ruled out. This conclusion might be 
avoided if the Z scale also decreases appropriately, we note. 

A second option to consider is the 2D equilibrium solution where, in the shorter-length 
case of zero Ai for example, |P| = r is a constant, say r — ► r 0 where from (4.3) rl = — Ci r /c 6r 
is positive and consistent with the value — cu/c 6i by virtue of the absorbed phase shift a. 
This 2D equilibrium is nominally an exact solution for all X, with zero induced-vortex 
flow, but it is unstable to 3D perturbations. For, if r = r 0 + 6r\, 0 = 861 with 6 small, then 


(4.1), (4.3) yield the linearized system 

Krr - 7W l!c = 2r„r,j/V 2 , (4.6a) 

r iFF“ Yr ik — ~ w i ?■> (4-66) 

t \y °, w i ~ 2r o r iz{-^F + <t>) as F -> oo, (4.6c) 

0 = Cl r T\ + C 3 r A 3 ifo + C 4 r* \ZZ izz ”b 3Cs r rgri, (4.6 d) 

0 = Ci,Ti + c 3 iA 3 iro + c^r^Yz ~ c ir^izz + 3cs«ro r i> (4.6c) 

with {W,t) * 6{W u n), ^31 = ri(X,0,Z), and Wi,ti to vanish as Y — »■ oo. The solution of 
(4.6a-c) may be decomposed in the form 

ri = ri exp (i(3Z + 7 X ) (4.7) 

and similarly for W\, Ti, A 31 , 0i and then the vortex solution is given by 

W-i — 2i/3r 0 riM(y),y = 7 *Y, 
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(4.8a) 


M{y) = B x Ai(y) - iny - Ai(y) Ai 2 {y x ){J^ y 2 Ai (y 2 ) tny 2 dy 2 } dyi 

h = r or i N(y), 

N(y) = B 2 Ai{y) + M'(y) + y' 1 + Ai{y) At~ 2 (yi){ J Ai(y 2 )dy 2 }dy x (4.86) 

where B x Ai( 0) = fan 7 /)/ 3 + <f> and B 2 Ai'( 0) = l/(3At(0)). Hence A 31 = 2y?^ 2 r 0 r 1 iV(0) 
with N( 0) = — e x — e 2 in*j and e x = 1.21 + 0.73 <f>,e 2 = 0.243, after some working. Then 
(4.6d,e) yield two coupled equations for r x ,0 x , and for these to have a nontrivial solution 
the eigenrelation 

/? 2 [( c 4 r + c ri) + 2 r^(c 3 r C 4 r + C3.C4i)T*(«i + 62 ^ 7 )] = 2ro(c 5 r c 4r + c 5 t c 4 .) (4.9a) 

must hold, determining the 3D spatial growth rate 7 in terms of the input 2D amplitude 
r 0 and the scaled spanwise wavenumber With the numerical values in (4.2) inserted 
into (4.9a) we then have the relation 

T*(« 1 + e 2 ln*f) = —d x /r\ - d 2 //? 2 (4.96) 

for 7(ro,/?) where the constants e x ,e 2 ,d x ,d 2 are all real and positive. See Fig. 2. Positive 
growth factors 7 exist for all nonzero ro, /?, the maximum growth factor being for large 
r 0 ,/? where 7 approaches the value exp(—e x /e 2 ) from below. Thus secondary 3D instability 
occurs for the 2D equilibrium state. 

The secondary 3D modes given by (4.9b) are completely absent in the undisturbed state 
of steady Blasius flow, we observe, since then (4.9a) is simply replaced by /? 2 (c 2 r + c 2 ,) = 
( c ir c 4r + c n c 4») or ?i = 0. The above secondary instability and those to be found in 
Sections 5,6 form some connection with Squire-mode destabilization (e.g., Herbert 1984, 
and see below), and the present one applies similarly to the nonparallel case Ai ^ 0 and 
rules out the 2D possibility as a terminal response of the full interaction. In addition the 
dependence of on /? tends to suggest that some spanwise focussing may be significant in 
the full nonlinear system (4.1). 

Third, there is the possibility of the nonlinear interaction continuing to downstream 
infinity in an algebraic form. The orders of magnitude present then suggest that 

W ~ X~*W + • • • , (r, A 3 ) ~ X{t , A 3 ) H , (4.10a, 6) 

r = |P| ~ X-»(enX)-&(z),0 ~ B(z) + • • • (4.10c, d) 

as X — ► oo, with the Z scale contracting such that Z = X~*z with z of 0(1). Hence the 
vortex equations become, with tj = Y / X* of order unity, 

d 2 W/dr 7 2 + ( y 2 /3)dW/d V - {rjz/2)dW /d~z + fa/6 )W = 0, (4.11a) 
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d 2 r/dr) 2 + (rj 2 /3)dr/dr} — (r7z/2)dr/dz — tjt = —dW/dz, (4.116) 

subject to 

W = [d(f 2 )/dz]/3,df/dr] = 0, r = A 3 at r? = 0, (4.11c) 

W — ► 0, r — ► 0 as rj — * oo, (4. lid) 

and coupled with the TS equation 

0 = (C2A1 + c 3 A 3 )P + c 4 d 2 P/ dz 2 (4.11c) 


where P = f exp(t'0). The TS-forcing of the vortex motion now appears only through the 
inner restraint on W in (4.11c), while the vortex forcing on the TS wave is due to the A 3 
term in (4.11e), and it is seen also that the amplitude-cubed TS effects become secondary in 
this type of spatial evolution. The solution of the vortex part (4.11a-d) leads to an induced 
skin-friction factor A 3 which enables (4. lie) to be written as an integro-differential equation 
for P(z), 


0 = 


c 2 A x — c 3 7z 



pm* 


dt 


P + C 4 P" 


(4.12a) 


~P=\P\ 2 , (4.126) 

where 7 = 2/[3tr(l/3)] is a positive constant and symmetry about z = 0 is assumed for 
convenience. The next question is whether an acceptable solution exists or not. Finite- 
difference computations have been performed for (4.12a,b) and these all suggest that there 
is no solution smooth for all z. Analytically, the large — \z\ behavior poses difficulties 
since the main suggested asymptote there, p ~ k s z 2 with fc 3 a positive constant, makes 
the integral in (4.12a) positive, so that the square-bracketed term cannot tend to zero 
then. Instead the solution appears to hit a singularity at finite z = z 0 — in which p tends 
to infinity like (z 0 — z)~s, and this seems in line with the computational findings. The 
singularity cannot be smoothed out locally, it seems, and so option 3 fails in that sense. 
The fourth option is that the nonlinear interaction (4.1) terminates in an algebraic 

A A A 

breakdown at a finite distance X , say as X — ► X $ — . The vortex part suggests that, in 
terms of order of magnitude, Y ~ ( X t — X)*,W ~ r 2 ~Z~ l ,t ~ A 3 ~ r 2 Y 2 Z ~ 2 and hence 
the dominant effects in the TS part are due to the c 3 , c 4 , C5 contributions, giving terms 
~ r 3 (X a — X)*Z *,rZ 2 ,r 3 respectively. For these to be comparable Z ~ (X t — X)s ~ r -1 , 

A A 

This indicates that the breakdown as X -» X t — takes the form 


W~{X.- X)- l L~*W , (r, A s ) ~ (X s - X)~H~\t, A 3 ), 
r = |P| ~ [X, - X)-*L-b[z),0 ~ 6(z ), 


I 

i 
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(4.13a, 6) 
(4.13c, d) 


to leading order, focussed near Z = Z, say with 

7-7.= {X t -X)*Lh, (4.13e) 

and Y = (X, — X)ai),L = — Ln{X t — X), the extra logarithmic factors here and in (4.10) 
being required by the inner condition in (4.1c). The coupled local governing equations and 


constraints are, from (4.1), 

d 2 W/dfi 2 - (fi 2 /3)dW/df) - ( fiz/3)dW/dz - f}W = 0, (4.14a) 

d 2 f/dfi 2 - (r) 2 /3)af/dr) - ( fiz/3)dr/dz - (2r)/3)f = - dW/dz , (4.146) 

W (oo, z) = f(oo ,z) = (df /dfi)(0, z) = W (0, z) — (r 2 )' / 3 = f (0, z) — A 3 = 0 (4.14c) 

for the vortex motion, and 

0 = c 3r A 3 r + c 4r (f" - rO n ) - C4i(2r'0' + rd") + c 5r f 3 , (4.14d) 

0 = cs,A s f + c 4 ,(r" - rO n ) + c 4r (2r'£' + rd") + c 5 ,r 3 (4.14e) 


for the TS part. Here once again the TS-forcing on the vortex motion appears only 
through the inner condition, in (4.14c), with the vortex forcing on the TS wave appearing 

A 

via A 3 in (4.14d,e). In contrast with the previous possibility however the amplitude-cubed 
TS contribution here is a primary effect. The solution of (4.14a-c) can now be used to 
express A 3 as an integral involving (f 2 )" and from this and (4.14d,e) we obtain the integro- 
differential equation 

0 = c 3 f )iz~ 1 J ^71 } P + c 4 p" + c 5 p|p| 2 (4.15) 

with p = |p | 2 and p(z) — rexp(i0), while 7 = 3?/r(l/3) is a positive constant and again 
symmetry about z — 0 is assumed. The issue therefore depends on whether (4.15) admits 
an acceptable solution or not. We note that as \z\ — *■ 00 matching could be achieved with 
the flow solution outside, since then (4.15) suggests the behavior p oc \z\~ l and hence 
P, W, etc., become 0(1) when Z — Z t is 0(1). Finite-difference computations performed 
for (4.15) suggest the same conclusion as for option 3, however, namely the occurrence 
of a singularity at finite distance, this time with p tending to infinity like (zo — z)~% as 
z — > Zq — , say. Again the singularity seems non-removable, thus ruling out option 4. 

The third and fourth possibilities above both represent the formation of concentrated 
tongues of high vorticity and/or TS amplitude. 

Option 5 for the present interaction is that the zero state is approached far downstream, 
with 

{W, t, P) 0 as X -+ 00. (4.16) 
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The stability of this state can be analyzed as for option 2, thus leading to (4.6) in the 
zero - Ax case but with ro replaced by zero. Hence at large X say A31 is zero as well and 
no growth of the perturbation is possible as anticipated in earlier comments. The option 
(4.16) means that the Type I interaction peters out but that may in turn allow the Type 
III interaction addressed in Section 6 below to be set off subsequently further downstream, 
since the Type III involves smaller TS amplitudes. 

The sixth option to consider has a finite-distance singularity in which the TS pressure 
remains finite but exhibits a singular behaviour in its gradient, as distinct from the strong 
singularity proposed in option 4. This option follows from a simple exact solution which 
is a generalization of the 2D one mentioned earlier to include a single spanwise-mode 
dependence, P = P 0 (X) exp(t/?Z) or r = R o (X),0 = 0o(X) + PZ, say. Here there is no 
induced vortex W — t flow, only the free W — r flow due to the input upstream, giving the 
skin-friction effect A 3 (X,Z) in general. The solution fits together if A 3 is independent of 
Z since then Ro satisfies 

0 = [cir — C4 rP 2 + C2 r AiX + C3 r A 3 (X)j i?o + C$ r Ro. ( 4 - 17 ) 

A ^ 

Hence as X increases the amplitude Ro can hit zero in a square-root fashion if A 3 (X) goes 
sufficiently positive that the square-bracketed term above approaches zero. An extension 
to the general case of (4.1) may be made next (see also Fig. 3). The proposal then is that 

A A A 

as X — * X,— for some finite station X s the irregular response 

P = P 0 (Z) + {X, - X)^P l (Z) + ( X t - X)P 2 (Z) + • • • , (4.18a) 

A 3 = A 30 (z) + {X, - X)X 31 (Z) + 0(X, - X) ? (4.186) 

is encountered. This is consistent with (4.1) provided that the successive equations 

0 = [c x + c^X-iX, + c 3 A 3 o(^)]-Po "b c *Pq ~b c 5 P 0 \Po\ 2 (4.19a) 

0 = (cx + c 2 XiX, + c $ X i0 {J)}P 1 + c<P” + C5[2 |Po| 2 Pi + PqPi] (4.196) 

0 = [cx + c 2 XiX, + c 3 A 3 o(Z)]P2 ( — ^ 2 Ax + c 3 A 3 i(Z)]Po 

+c 4 P 2 w + c 6 [2|P 0 | 2 P 2 + P 0 2 P; + Px 2 P 0 * + 2P 0 |Pi| 2 ] (4.19c) 

hold for Po,PijP 2 in turn. Here the vortex skin-friction factors A 30 , A 31 remain arbitrary 
in the sense that they depend on the history of the flow upstream, and this arbitrariness 

makes the solution of (4.19a-c) difficult in general. As a model however suppose that A 30 

is constant and the coefficients c n (n = 1 — ► 5) are real. Then a solution of (4.19a) has 
P 0 equal to a real constant, while (4.19b) allows P 4 cx sin(/fe) with P 2 c 4 = 2 c 5 P 2 and a 
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solution with the same periodicity requirement in Z can be found for P 2 in (4.19c). Hence 
the terminal form then is self-consistent. This model can also be enlarged to allow the c„ 
values to be complex as in (4.2) , and other models support the proposal above. Another 
enlargement possible has spanwise focussing where the Z-scale shrinks like a power of 
X s — X at the singular station and this appears to fit together as well. 

The overall interpretation for this option 6 is that the induced-vortex skin friction 
can act to produce an irregular response in the TS amplitude, without a corresponding 
irregularity in the vortex motion. Further, the dependence on the history of the induced- 
vortex motion emphasizes the role of the upstream starting conditions in determining 
whether option 6 arises or not. 

The computations referred to earlier seemed to confirm the availability of option 5, 
but they were not conclusive with regard to option 6. In particular it proved difficult to 
distinguish between numerical divergence of the iterative procedure used and appearance 
of a genuine square-root irregularity in the solutions. Failing other options (1-4 above 
and certain other ones tested), we tend to the view however that options 5,6 are the only 
ones attainable via the nonlinear interaction (4.1). The other options are associated with 
by-pass transitions: see Section 7. 

The spanwise dependence in option 6 again allows streak-like formations to appear, 
although less firm than in the options 3,4 in Hall and Smith’s (1988a) interactions and in 
the Type II interaction below. Option 6, with its shortening length scale streamwise, can 
lead into the type II interaction, which is studied next. 

5. “TYPE IP INTERACTIONS 

The earlier arguments for Type I interactions may be modified to account for an alter- 
native, Type II, interaction as follows (again see Fig. 1). The relative errors referred to at 


the start of Section 2 have the orders 

/3 2 [from the warping], (5.1a) 

h 2 [from amplitude-cubed feedback], (5.16) 

/ 3 2 h 2 k~ * [from the induced vortex shear], (5.1c) 

h 2 [possible from streamwise TS modulation], (5. Id) 

k [from the streamwise vortex modulation], (5.1e) 


and in the Type I case above the balance of (5.1a-c) holds, yielding (2.2). The present 
alternative interaction (Type II) occurs for /3 ~ h again, balancing the effects (5.1a,b), 
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but k ~ h? now, so that the typical streamwise spatial evolution is faster here. Then the 
vortex-shear effect (5.1c) is of order h% and is relatively small, giving a reduced vortex 
influence on the nonlinear TS evolution. 

The Type II interaction therefore has 

/ 3 ~ h,Jc ~ h?,h << 1 , (^* 2 ) 


and it can be seen from modification of the working in sections 2,3 that the new sizes (5.2) 
lead to a TS amplitude equation in which dP n /dXi replaces the interaction term oc A 3 
on the left-hand side in (2.14) and absorbs the er-term there also. Hence we are left with 
solving the nonlinear equation 


— ( c i + ci\iX)P + C4 2 + c 5 P\P\ 2 (5.3) 

oX dZ 

instead of (4.1e), the vortex effect (4.1a-d,f) then being of secondary significance. The 
constants Cj,C2,C4,C5 are as in Section 4. 

The governing equation (5.3) for the TS pressure P(X,~Z) is a generalized cubic 
Schrodinger or Ginzburg-Landau equation, and its main properties seem to be clear. 
First, the pure 2D version where the C4 term is suppressed produces a supercritical bi- 
furcation, namely that in Smith (1979), Hall and Smith (1984), leading to a saturation 
amplitude downstream as X increases: see below. Second, however, the 3D version with 
Z-dependence present shows that there is unbounded secondary instability of the 2D state. 
This is due mainly to the coefficient c 4r begin negative. Thus in polar form P = rexp(iO), 
which replaces (5.3) by the two real equations 


r x ( Clr c 2 rAiX)r + c ir (r zz r0j) z + rO zz ) + cg r r 3 , (5.4a) 

r0% = (ci,- + c 2i AxX)r + c 4 ,(r I j - rOj) + c^rjOj + rQ z z ) + c 5i r 3 , (5.46) 
the pure 2D flow solution has (r, 6 ) = (r 0 , #o) (-AT) where 

r 0 — (cir d - c 2 rAi-AT)r 0 d - c$ r r q, and Qq — (ck -\- C 2 ,Ai^) -|- C 5 ,rQ, (5.5a, 6) 


so that with k = c 2r Ai positive (the nonparallel-flow case) 


r 0 = 


2 [X _ 

exp(/cX )/{— 2c 2r /_ exp(Kq 2 )dq} 
J Xq 


, X = X + Cj r /c 1 , 


(5.5c) 


for X beyond an initial station downstream of Xo, while if Ai is zero (the parallel-flow 
case) 

r = [c lr ai exp(2c lr X)/{l - aiC 5r exp(2ci r X)}]£ (5.5d) 
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where the constants <ii,Ci r are both positive. The saturated 2D solution downstream 


exhibits parabolic growth ro ~ (/cX/(— C 5 r ))* or tends to the constant value r 0 = = 

(cir/(— c 5 r))^ as X — * oo. A small 3D perturbation of the 2D solution, 

M] = [r 0 A](X) +6{r 1 ,O 1 }{X)exp{i0Z) + 0(6 2 ), (5.6) 

with spanwise wavenumber /?, is then controlled by the linearized system 

r[ = (ci r + c 2r \iX - ~jf c 4r )r t + c 4 ,r o 0i + 3 c hr rlr u (5.7a) 

r 0 6[ = ~p 2 c 4i ri - p 2 c 4r r o 0i + 2c 5i rlr u (5.76) 


from (5.4a, b). As posed, the secondary instability problem (5.7a,b) represents a nonparallel 
marching problem in X for the spatial evolution of the 3D perturbation, from given starting 
conditions at some finite X and with r 0 {X) specified in (5.5c or d). Downstream at large 
positive X however the 2D equilibrium r 0 — > Too (constant) holds if we focus on the case 
of zero Ai and there the solution of (5.7) is sought in the form 

(ri,0i) = (f!,£i)exp(QX) (5.8) 

with the spatial-growth factor Q to be found. In effect, each of the primes in (5.7a,d) is then 

A A A 

replaced by a factor Q,r „> replaces ro(AT), and so substituting for 0\ from the equivalent 

A 

of (5.7b) into (5.7a) we obtain, for nontrivial solutions, an eigenrelation determining Q in 
terms of /7, which gives 

Q = c 5 rr^ - c 4r f ± [ c\ r r + 2c 4 ,c 5 ,r^^ 2 - c 4 ,^]?. (5.9) 

The maximum spatial growth rate Q r therefore arises at large spanwise wavenumbers /?, 
for which Q r ~ — c 4r ^ is large and positive. Indeed, the same large unbounded growth 
also occurs at finite X, within the framework of (5.7) (with Ai zero or nonzero), as a 
pronounced short-wave instability. This again raises the possibility of spanwise focussing 
taking place in the nonlinear system (5.3). 

Finally, the fully nonlinear 3D version, in (5.3), produces in general a break-up of the 
nonlinear TS solution within a finite distance, associated with the formation of “vorticity 

A A 

tongues.” The break-up, at the station X — ► X,— say, has 

\P\=r~{&,-&)-h{Z) + — , 

arg P - 0 ~ kitn(X s - X) + 0{Z) 

Z-Z t = (X t -X)^Z, 
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(5.10a) 

(5.106) 

(5.10c) 


focussed at a particular spanwise station Z = Z s : see Fig. 4. Here ki is a constant. So 
the governing equations (5.3) or (5.4) reduce to the nonlinear ordinary-differential system 


i(f + Zr') = c 4r (f" - rb n ) - c 4 ,(2f'0' + rb ") + c 5r f 3 , (5.11a) 

r(—k\ + \Z0') = c 4 ,(f" - rb 12 ) + c 4r (2r'0' + rb") + c 5 ,f 3 (5.116) 

It 

for f, b. The existence of acceptable solutions of this terminal form (see also Hocking and 
Stewartson (1972) and references therein) is suggested simply by a linearized version for 
example in which the cubic terms in (5.11) are negligible. One solution then has 


f oc exp[-£ 2 (-c 4r )/4|c 4 | 2 ], 
b = constant — Z 2 c 4l /4|c 4 | 2 , 


(5.12a) 


and the other solution is expressible in terms of error functions, giving an algebraic 0(Z~ 1 ) 
decay in f at large \Z\ which matches with the flow solution outside where r is 0(1) at 0(1) 
values of Z — Z t . Both of the above solutions are acceptable because the coefficient c 4r is 
negative. 

The finite-distance break-up in (5.10) - (5.11) is clearly a strong one and it may be 
interpreted as producing strong tongues of concentrated streamwise vorticity and TS am- 
plitudes in the singular form (5.10). New higher-amplitude shorter-scale phenomena must 

A ■■■■ 

then come into operation, closer to the break-up point X t ,Z t , to continue the development 
of the tongue. Specifically, the full 3D nonlinear system (2.1) is then reinstated. This is 
mainly because the pressure amplitude p due to (5.10) is becoming of order h(X , — X)"*, 
which reduces 0(1) when the streamwise distance \X — h~ 2 X s \ reduces to 0(1), from (5.2). 
Simultaneously the Z-scale contracts to 0(1), from (5.1) with (5.10), the induced-vortex 
strength rises to 0(1) also, and of course an 0(1) time scale T is present because of the 
primary wave (2.9). Hence the further evolution of this tongue or streak is controlled by 
the fully nonlinear triple-deck system (2.1). Some computational solutions for such 3D 
unsteady triple-deck flows are given by Smith (1988a), while the likelihood of finite-time 
break-ups occurring even in the full triple-deck framework is shown by Smith (1988b). 


6. “TYPE IIP INTERACTIONS 

Here once again the reasoning concerning the error sizes in (5.1a-e) may be modi- 
fied/extended to describe a new kind of nonlinear interaction possible, Type III (see also 
Fig. 1). This last type has j3 being of order unity, i.e. faster spanwise dependence than 
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before, but with h still being small and k again being 0 (h 3 ) as in the first type of interac- 
tion. So now the vortex-shear influence (5.1c) is also of order unity, balancing the warping 
effect (5.1a). 

This Type-Ill nonlinear interaction, in which 

/?~l,fc~/i 3 ,/i«l (6.1) 


therefore provokes a nonlinear effect on the whole flow, in particular completely altering 
the mean-flow profile., despite the smallness [h « 1] of the input TS wave. In more detail, 
some allowance has to be made for the logarithmic interplay (see Section 2) between the 
buffer and the near-wall zones, and in fact the near-wall zone has the expansion 


j — 


( u — A„Y + h'{uiE + c.c.) + h l2 u v + 
v = h'(viE + c.c.) + h l2 v v + • • • , 
w — h'{wiE + c.c.) + h l2 w v + • • • 


(—tnh) * 

[ w = h'(wiE + c.c.) + 
now, while in the buffer zone where Y = h~ x Y we have 


(6.2a) 

(6.26) 

(6.2c) 


u = h~ 1 U v (X,Y,Z) + •••, (6.3a) 

v = hV v (X,Y,Z) + .*•, (6.36) 

w = h 2 W v (X,Y,Z) + •••. (6.3c) 


In (6.2a), X V (X,Z) is the unknown vortex-induced shear, and a short- and a long-scale 
streamwise dependence hold, in that 

d d 3 d 

ax ax + h dx‘ 

The subscript v refers to induced-vortex quantities. Also the unknown pressure and dis- 
placement take the forms 

p = h'{PE + c.c.) + •••, (6.5a) 

A = h- 1 A v + h'{AE + c.c.) + ••• (6.56) 



where, as in (6.2), the smallness of the input TS disturbance is evident. 

The principal feature of the near-wall-zone form in (6.2) is that it involves a small 
perturbation of the unknown shear flow u = A V Y and as such it leads at first order to the 
linear equation 


d 2 P 1 d\ v dP 
dZ 2 X v dZ dZ 


J — ol 2 P = AQ 


(6.6a) 
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for the TS pressure P, after some working as in Smith (1980). In (6.6a), A may be replaced 
by a linear integral of P, stemming from the interaction law (2. If), or the equivalent upper- 
deck-flow setting may be invoked, i.e. solve 


(< dl + d*-a 2 )p=0 , (6.66) 

p = P, d^p = -a 2 A at y = 0+, (6.6c) 

with p bounded at infinity, to obtain the P — A law. Again, with some notation of Section 

2 , 

o t a;i / c \ a jt 

(6.6 d) 


j _ £ , £o-At 

“2 2 Ai 


6>At{, / £qk\ 


and a, fl must be real. We observe that (6.6a) looks linear but strong nonlinear dependence 
is still present, nevertheless, through the unknown shear A„. At second order a vortex-flow 
contribution (among others) is induced due to amplitude-squared forcing and that yields 
the logarithmic growth 


w v ~ -\- 2 d z [\P\ 2 + a- 2 P z P* z ]tnY as Y -* oo 


( 6 . 7 ) 


similar to the growth described in Section 2, with the subscript v again referring to the 
vortex component. 

The growth (6.7) and the corresponding growths proportional to Y z tnY, Y InY in u v , v v , 
respectively, are responsible for the considerable vortex effect set up in the buffer zone as 
described by (6.3). The controlling equations in this buffer are the full nonlinear 3D 
boundary-layer equations. 


du v dV v dw v _ 
dx v dY ^ dz 

(6.8a) 

u v du v v v du v w v du v _ d 2 u v 
dx h d¥ + dZ d ¥ 25 

(6.86) 

u v dw v v v dw v w v dw v _ d 2 w v 

dX dY dZ dY 2 

(6.8c) 


from substitution of (6.3) into (2.1a-c). It is noteworthy that the induced pressure gradient 
for the vortex motion is negligible here. This is because the main displacement (in (6.5b)) 
associated with (6.3), although of order h ~ l , depends mostly on the slow variable X, not 
the fast one X, and hence the corresponding induced pressure p is only 0 (h 5 ), from (2. If) 
or from analysis of the upper-deck behavior of the vortex flow, which then provokes a neg- 
ligible feedback in the momentum balances (6.8b, c). Thus the outer boundary conditions 
here are 

U V ~Y + A V (X , Z), W v -+ 0 as Y -»• oo (6.8d) 


22 



where the displacement function A V {X, Z) is to be found from the solution of (6.8a-c). 
The inner boundary conditions account for the match with (6.2) combined with (6.7), so 
that 

U v = V„ = 0 ,W V = -\~ 2 d z \\P\ 2 + cT 2 P z Pz ] at F = 0. (6.8e) 

Finally here, the unknown vortex shear 

A„ = ^(X,0 ,Z) (6.9) 

links the buffer and near-wall solutions together, providing the overall nonlinear interac- 
tion between (6.6), (6.8). Another noteworthy point, to repeat, is the smallness of the 
TS amplitude compared with the considerable vortex motion that it provokes; and the 
TS-squared forcing of the vortex motion in the buffer zone is dominated by the effec- 
tive spanwise-slip- velocity condition in (6.8e); other forcings such as those through the 
momentum equations are of less importance, at least at the current stage. 

The full Type-Ill interaction is controlled then by (6.8a-e) subject to (6.6), (6.9), for 
the driven vortex flow and the driven TS pressure, respectively. 

Secondary instability of a 2D TS input can be established as for the previous Types I, II 
and as for the channel-flow interaction of Hall and Smith (1988b). The secondary instability 

3 

of the pure 2D state, i.e. P = Po(X) with U v = Y, V v = W v = 0 and a = a 0 = dj , A„ = 1 
(which gives an exact solution of the system), arises as a small 3D perturbation in the 
form 

[P,A,U v ,V v ,W u \ = [P b,A 0 ,F,0,0] 

+6[PW cos PZ, A™ cos PZ, UW cos pZ, V™ cos PZ, sin PZ] + • • • (6.10) 

and A„ = 1 + ^A^ 1 ) cos PZ + • • •, with 8 being small. Here we take the case of Po,A 0 being 
constant, within the present length scales, and all the Z-dependence is as shown explicitly 
in (6.10). The 3D induced- vortex flow is therefore controlled by the linearized equations, 
from (6.8a-e), 


ug + f£> + pwW = 0, 

{6.11a) 

Yu^ + V™ = U%, 

(6.116) 

r«j? = w vf’ 

(6.11c) 

subject to the constraints 


tfW 0,IFW -» 0 as F-+ oc, 

(6.11<i) 

U {1) = F^ = 0, = p[P 0 P w * +' P 0 *P (1) ], 

(6.11e) 
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with the last condition representing the TS-forcing effect. Along with this, the TS pertur- 
bation satisfies, from (6.6), 

(P* + a 2 0 )pM = a *AW + XX^Pocto 1 , ( 6 . 11 /) 

(<*o + /3 2 ) *P (1) = a£A (1) . (6.11gr) 

Here (6.11g) follows from (6.6b, c) and P 0 = a 0 A 0 , while X = [| — 2 ^*° (1 + i^)] a o = -a o^ 
£o, c are defined just prior to (2.12b) and in (2.15), respectively, and the vortex-forcing term 
is A^ = dUW ( X , 0) / dY . The whole 3D perturbation is governed by the interactive system 
(6.11a-g) or by its integro-differential counterpart analogous to that in Appendix A. So 
the spatial development of the perturbation with downstream distance X is in general a 
nonparallel one, starting from input conditions at some finite station upstream. We turn 
therefore to the far-downstream response where spatial instability is proposed, 

P W „ pM exp(sX) asX^oo, (6.12) 


say, and similarily for the other variables, with the constant growth factor s to be found. 
Then the vortex part (6.11a-e) yields 


WM = K (e) dl^l ,q = sSY, 
A i(0) 

duM _ -pK^Ai'in) 

W ~ a lA»(0) 


with = P[P 0 PM* + c.c.], so that the vortex skin-friction factor is 

A (e) = Z^M [PoP ^ + p*p(*)], 

ssAt(O) 


(6.13a) 

(6.136) 


(6.13c) 


Hence the TS part (6.11f,g) reduces to the complex linear equation 

{P 2 + a 2 )pM = a 0 {a 2 0 + p 2 ) *P« + [ P 2 P (e) * + |Po| 2 P ( Ol (6.14) 

a 0 s® J 

where T = — At'(0)/At(0) is a positive constant. The real and imaginary parts of (6.14) 
then provide an eigenrelation for the determination of s, giving the result 


s* = 2Ta 0 1 )l r \Po\ 2 p 2 {a 2 + p 2 ) »[(o$ + $)» - a 0 ] (6.15) 


So the main properties of the growth factor s are that, first, it is real and positive, since 

)l r > 0; second, s decreases monotonically with increasing spanwise wavenumber /?, from 

2 

the 0(1) value s 0 at P = 0+ to |s 0 as P —► oo, where Sq = 4r^/ r |P 0 | 2 , and hence the 
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fastest growth is for nearly planar secondary modes, cf. Sections 4,5; third, s increases 
monotonically with |P 0 |, like | JP 0 | 3 ; and, fourth, this secondary growth is present if and only 
if the input 2D TS amplitude |P 0 | is present. Other aspects are similar to those in Section 
4 including the spanwise focussing associated with the large-/? instability for initial-value 
problems. 

Secondary instability of a similar kind also occurs at relatively high frequencies. 

In the nonlinear regime there are several options which could be put forward for the 
ultimate behavior of the full interaction (6.6), (6.8), (6.9) as the downstream distance X 
increases, some of the options resembling those in Section 4. Again, the Type-III interaction 
connects up with a global-scale nonlinear interaction currently being studied by the author 
and Professor P. Hall. The Type-III interaction involves short-scale/long-scale balancing 
(via (6.4)) and it is clearly a robust one in that it changes the mean-flow characteristics (see 
(6.3)), even though the TS amplitude is still small (see (6.2)). Its full nonlinear properties 
remain to be studied. 

7. FURTHER COMMENTS 

We finish by noting the following items (a) -(f) on the nonlinear TS/vortex interactions 
studied above. 

(a) Each of the three Types I- III of interaction addressed in this work can be triggered 
by the input disturbance upstream, depending on the latter’s amplitude spectrum as in- 
dicated by the sizes specified in (2.2), (5.2), (6.1) [see also (d) below]. Alternatively, each 
interaction can start as a form of secondary 3D stability of the 2D TS input upstream. 
Beyond that however nonlinearity takes control. 

(b) There are some connections between the three Types I-II, as their governing equa- 
tions in Sections 4-6 suggest, and a match can be established also with the Hall and Smith 
(1988a) study concerning oblique TS waves. Further, the main options 5,6 for the behavior 
of the Type I nonlinear interaction may lead on into the Type III, II interactions, respec- 
tively, downstream. Type II has a faster streamwise response than Type I, while Type III 
has a faster spanwise dependence: see the scales in Sections 2,5,6. 

(c) The Types II, III appear to be much more powerful and dangerous nonlinear inter- 
actions than Type I. Type III involves only a small TS amplitude (see (6.2)) but despite 
that the mean-flow profile is completely altered (see (6.3)), due to the induced vortex 
motion. Similarly, Type II starts with small TS amplitudes (see (5.2)) but the ensuing 
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finite-distance break-up then causes the full triple-deck system (2.1) to be triggered [Smith 
1988a, b], thus completely altering the mean-flow profile again. 

(d) By-pass transition processes are also possible and of much interest. They are 
associated with other scales of input upstream, say, and in principle these can activate 
all the options 1-6 described for the Type-I interaction and the corresponding options 
available for Type III as well. These and other subsequent stages seem to merit further 
research. 

(e) The extra effects of wall curvature, cross flow and compressibility on the nonlinear 
interactions remain to be studied. All three can be captured by modifications of the triple- 
deck starting point (2.1) for instance: see earlier studies of compressibility effects in Smith 
(1987), cross- flow effects in Stewart and Smith (1987) and wall-curvature effects in Hall and 
Smith (1987, 1988b), among others. In turn, any of the three effects could be incorporated 
initially as an extra contribution to the nonlinear interaction equations of Sections 4-6, 
and it would be interesting to see their influence on the interactive flow properties. 

(f) No quantitative comparisons with transition experiments [see earlier references] have 
been attempted yet, but the possibility of a reasonably firm link exists there in view of 
the predicted formation of streak-like behavior (among other phenomena), with spanwise 
concentrations of the vortex flow and the TS amplitudes (Sections 4-6). 
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APPENDIX A: AN INTEGRAL FORMULATION FOR THE TYPE I IN- 
TERACTION 

An integral formula for the vortex skin-friction factor X 3 (X,Z) can be obtained from 
the vortex-flow equations (4.1a-d). A Fourier or Laplace transform in X, for instance, 
enables the solution for the transforms of W, r to be expressed in terms of the functions 
M,N introduced in (4.8a, b). Hence we find the result 

/* liftf. Z)(*- er f [Si + g,ln(X - ()]d( (Al) 

aZ J *o 

in the case of a vortex flow starting from rest, or zero shear, at a station X = X 0 , say. 
Here the constants gi,g 3 are given by 

a . = -«»/r Q) ,n = { e i - sir- (|) } /r(|) (A2) 

with ei,e 2 specified just after (4.8b). 

So the nonlinear Type-I interaction can be written concisely in an integro-differential 
form, by coupling (Al) with (4.1e). This reproduces the governing equations for all the 
options 1-6 presented in Section 4. A similar formulation can be constructed for the 
oblique-wave/ vortex interaction in Hall and Smith (1988a) and is used in the wave/ vortex 
interactions being studied by Mr. N. D. Blackaby and Mr. P. A. Stewart. 
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FIGURE CAPTIONS 


Fig. 1. 

Fig. 2. 
Fig. 3. 
Fig. 4. 


Nonlinear interaction Types I-III between long-scale vortices and short- 
scale waves, in normalized coordinates, for waves of low input amplitude 
0(h). The values of the length scales /3~ 1 ,k~ 1 for I-III are given in the 
text (sections 2,5,6 in turn). 

The scaled 3D spatial growth rate 7 versus spanwise wavenumber (3 for 
various input 2D amplitudes r 0 , in secondary instability for Type-I inter- 
action. 

The square-root spatial breakdown of the wave amplitude (option 6) for 
Type-I nonlinear interactions, due to the vortex-shear effect (A 3 ), leading 
to shortening of the interaction length scale. 

Nonlinear breakdown of the Type-II interaction, via spanwise focussing 
and spatial blow-up of the wave amplitude. 
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Spanwise wavenumber, (3 
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